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ABSTRACT 

V471 Tauri, an eclipsing system consisting of a hot DA white dwarf (WD) and 
a dK2 companion in a 12.5-hour orbit, is the prototype of the pre-cataclysmic 
binaries. The late-type component is magnetically active, due to its being con- 
strained to rotate synchronously with the short orbital period. During a program 
of ultraviolet spectroscopy of V471 Tau, carried out with the Goddard High Res- 
olution Spectrograph (GHRS) onboard the Hubble Space Telescope, we serendip- 
itously detected two episodes in which transient absorptions in the Si III 1206 A 
resonance line appeared suddenly, on a timescale of <2 min. The observations 
were taken in a narrow spectral region around Lyman-a, and were all obtained 
near the two quadratures of the binary orbit, i.e., at maximum projected sepa- 
ration (~3.3-R0) of the WD and K star. 

We suggest that these transient features arise when coronal mass ejections 
(CME's) from the K2 dwarf pass across the line of sight to the WD. Estimates 
of the velocities, densities, and masses of the events in V471 Tau are generally 
consistent with the properties of solar CME's. Given our detection of 2 events 
during 6.8 hr of GHRS observing, along with a consideration of the restricted 
range of latitudes and longitudes on the K star's surface that can give rise to 
trajectories passing in front of the WD as seen from Earth, we estimate that the 
active V471 Tau dK star emits some 100-500 CME's per day, as compared to 
~l-3 per day for the Sun. The K dwarf's mass-loss rate associated with CME's 
is at least (5-25) x 10~^^ Mq yr~^, but it may well be orders of magnitude higher 
if most of the silicon is in ionization states other than Si III. 
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Subject headings: coronae: stellar — mass loss — stars: V471 Tau 

1. Introduction 

It has long been difficult to estimate with confidence the mass-loss rates from cool dwarf 
stars. On the basis of X-ray emission, there is clear evidence that coronae exist around many 
cool dwarfs. Based on analogy with the Sun, it is likely that thermally driven winds emerge 
from these coronae. However, the column densities of these winds are so small that direct 
detection is usually very difficult. 

The thermally driven wind is, however, not the only means by which the Sun loses mass. 
There are also episodic ejections of material into the wind. These coronal mass ejections 
(CME's) contribute on average some 16% to the solar-wind mass flux at solar maximum 
(Jackson & Howard 1993). The CME's occur when magnetic configurations in the corona 
can no longer find an equilibrium with the ambient plasma. 

Because CME's are intrinsically magnetic in nature, it is expected that in stars where 
magnetic activity levels exceed those in the Sun, CME's should be larger in scale, and 
occur more frequently, than in the Sun. In particular, lower-main-sequence stars belonging 
to the category of flare stars are sites of particularly high magnetic activity. The energy 
released in flares on such stars may exceed by several orders of magnitude the energy that 
is released in the largest solar flares. In such magnetically active stars, we anticipate that 
CME's might contribute significantly more than 16% to the total mass flux. In this paper, 
we report direct spectroscopic evidence that points strongly to the existence of CME's from 
a well-known magnetically active star. 

Our target is the late-type flare star in the V471 Tauri binary system, a member of the 
Hyades star cluster. V471 Tau is an eclipsing binary with a period of 0.52 day, consisting of a 
cool main-sequence star (dK2) and a hot white dwarf (DAI. 5). Orbiting at only ~3.4 stellar 
radii from the dK star, the hot white dwarf (WD) provides a superb UV background source 
for probing the environment above the cool star's surface. In 1986, X-rays from the WD were 
found to be rotationally modulated with a period of 9.25 minutes (Jensen et al. 1986), and 
it was hypothesized that the modulation might be due to accretion of heavy elements onto a 
magnetic polar cap on the WD. Direct spectroscopic conflrmation of the accretion hypothesis 
was obtained when Sion et al. (1998) discovered a rotationally modulated, Zeeman-split UV 
metallic absorption line in the photospheric spectrum of the WD. 

The material that accretes onto the WD in V471 Tau is presumed to originate in a 
wind from the dK2 companion. Reliable evidence that the K2 dwarf is magnetically active 
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is provided by optical spectroscopic and photometric flares which are "common and intense" 
(Young 1976; Young, Rottler, & Skumanich 1991; see also Beavers et al. 1979), by X-ray 
flares (Young et al. 1983), and by radio flaring activity (Patterson, Caillault, & Skillman 
1993; Nicholls & Storey 1999). The high levels of magnetic activity in this particular K2 
dwarf are attributed to fast rotation: the star is tidally locked to the short orbital period of 
the binary. This leads to an angular velocity for the K2 dwarf, fldK, that is some 50 times 
greater than the solar value, O©. As a result, dynamo action is expected to be much more 
efficient than in the Sun. 

Associated with enhanced magnetic activity, wc expect an increase in the efficiency 
of mass ejections in the form of CME's. In support of this, wc note that Young (1976) 
has pointed out a possible connection between mass-loss events in V471 Tau and episodes of 
activity on the K2 star. Moreover, transient absorptions in Mg II in lUE spectra of V471 Tau 
were reported by MuUan et al. (1989), and interpreted as evidence for CME's. Based on the 
Mg II line strengths, MuUan et al. estimated the mass-loss rate associated with the CME's 
to be at least 10"^^ Mq yr'^ 

Eclipse timings show that the orbital period of V471 Tau changes with time (Tunca et 
al. 1979). If the period changes are secular in nature and are due simply to mass loss from one 
component, the rate of period change may be used to determine mean mass-loss rates. This 
argument led to an estimate of Med of order 10^^ Mq jt^^ for V471 Tau (Tunca et al. 1979). 
However, it is alternatively possible that the apparent period changes arc strictly periodic 
and arc due to light-travcl-timc effects in a triple system, as has been argued by Ibanoglu 
et al. (1994) and Guinan & Ribas (2001). Yet another possibility, perhaps the most likely, 
is that the orbital period modulations in V471 Tau and other active binaries arise from 
the magnetic activity. A specific mechanism is the gravitational interaction between the 
orbit and variations in shape of the active component in the system, as discussed by, e.g., 
Applegate (1992), who argues that it is plausible for the case of V471 Tau. In view of these 
alternatives, the cited value of Meci is probably a gross overestimate. In a later section, we 
will support this conclusion using an independent line of argument: we will show that the 
mass-loss rate from the K2 dwarf cannot exceed 10~® Mq yr~^. 

In the present paper, we report on UV observations of V471 Tau obtained with the 
Goddard High Resolution Spectrograph (GHRS) onboard the Hubble Space Telescope {HST). 
The spectra show clear evidence of transient absorptions in the 1206.51 A resonance line of 
Si III, which we attribute to CME's passing in front of the WD component. The Si III line 
is representative of hotter material than Mg II, and we therefore expect that the Si III line 
more nearly reflects the coronal conditions in V471 Tau. To the extent that this expectation 
is valid, the present study provides a more rehable analysis of CME properties than that 
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based on the earlier lUE data (MuUan et al. 1989). 

We note that V471 Tau is a detached binary, in which neither star fills its Roche lobe. 
Therefore, the mechanism that drives mass loss from the K2 dwarf has nothing to do with 
Roche-lobe overfiow. Instead, the mass loss from the dK2 star must depend for its origin on 
the same processes to which a single star has access. For that reason, we may meaningfully 
consider processes that are known to operate in the Sun in our attempt to interpret our 
observations of V471 Tau. 



2. Observations 

The GHRS observations that form the basis of the present paper have been described in 
detail by Sion et al. (1998) and O'Brien et al. (2001). The observing was optimized primarily 
to determine the radial velocities of the WD component, and also to search for spectroscopic 
variations on the rotational period of the WD; our discovery of transient CME events was 
thus entirely serendipitous. Briefly, in 1993, 1994, and 1995, in the course of 12 HST visits, 
128 individual spectra of V471 Tau were obtained centered on the wavelength of Lja. At 
these wavelengths, the hot WD completely dominates the spectrum, aside from a narrow 
chromospheric Lya emission line from the K dwarf. The dispersion was 0.07 A diode-\ 
and the wavelength range covered was 36 A, broad enough to include, by chance, the Si III 
1206 A resonance fine. In 1994-95, there were eight HST visits, during which each individual 
sub-exposure had a duration of 127 s, chosen in order to resolve the 9.25-minute rotational 
period of the WD into four time bins. During each visit, a series of 16 sub-exposures was 
taken, lasting (including overheads between sub-exposures) a total of 36.3 min. During the 
four visits in 1993, the sub-exposures were 101 s, and the total observation duration per visit 
was 29.0 min. Each visit was planned to occur near orbital quadratures, i.e., at the maximum 
projected separation between the WD and the dK2 components as seen from Earth. The 
total time that the star was under observation during the 12 visits was 6.8 hours. 

Transient absorptions in Si III A1206 A were seen twice during the course of the ob- 
servations. Fig. 1 shows a pictorial representation of the spectra obtained during the more 
dramatic event, which occurred on 1994 Oct 20. Sixteen individual 127-s exposures are 
stacked in Fig. 1 from top to bottom in order of increasing time, in order to simulate a 
single-trailed photographic spectrogram. Very broad photospheric Lya absorption wings 
extend across the spectrum, at constant strength in all 16 spectra, and there are emission 
features from the dK chromosphere and a deep interstellar Lyct absorption centered in the 
core of the hydrogen line. In the top four spectra, there is no sign of the Si III line, but 
this line appeared suddenly in absorption in the fifth spectrum. The absorption was strong 
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as soon as it appeared, and it remained comparably strong throughout the remainder of 
the 36-min series. This observation was taken at orbital phase 0.82 (WD receding and dK2 
approaching). When the next exposures at the same quadrature were obtained 3 days later, 
the Si III absorption had disappeared. As noted above, at this orbital phase, the projected 
separation of the stars is at its maximum of about 3.4 dK radii (O'Brien et al. 2001). 

A second, similar but less pronounced transient Si III event (not illustrated here) was 
detected near the end of another visit on 1994 Oct 17, at orbital phase 0.77. This transient 
made a sudden appearance in the 14th sub-exposure, and remained detectable during the 
final three sub-exposures of this visit. When the next visit at the same orbital phase occurred 
2 days later, the Si III feature had again disappeared. 

In Fig. 2, we plot the equivalent width (EW) of the Si III line versus time during the 
36 minutes of observation on 1994 Oct 20. The EW was below detectability (<40 mA) 
during the first four sub-exposures. Then, in the fifth sub-exposure, the line appeared at 
EW RillO mA. The line grew in strength to a maximum of ^275 mA in the eighth sub- 
exposure, and then fluctuated approximately between these limits for the remainder of the 
visit. 

The heliocentric radial velocity of the Si III line during the 1994 Oct 20 absorption event 
(determined by summing the 12 spectra showing the line, and measuring the wavelength of 
the hne center) was -1-176 ± 3 km s""^. At this time, the center-of-mass velocity of the WD 
was -1-188 ± 3 km s~^, according to the spectroscopic orbit of O'Brien et al. (2001). 

For the event on 1994 Oct 17 the Si III radial velocity was -|-220±10 km s^^, as compared 
with the WD velocity of +200 ± 3 km s^^ 

Thus both transient features had radial velocities within about ±20 km s~^ of the WD 
radial velocity. This indicates that, in the reference frame which rotates with the orbital 
motion, the transient absorber crossed the line of sight to the WD almost exactly transversely. 

We note that the sudden, strong absorption feature we discuss here is distinct from the 
much weaker and broader absorption Si III line arising in the photospheric surface spot on 
the WD discovered by Sion et al. (1998): the velocities are significantly different (note that a 
photospheric line would have an additional gravitational redshift of -1-49.9 km s~^), a surface 
"patch" of Si III ions would be rotationally modulated on the 9.25-min WD rotation period, 
and a photospheric line would show Zeeman splitting (see Sion et al.). 
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3. Interpretation of the Transient Events as CME's 

We will now discuss the transient absorption events in terms of CME's passing in front 
of the WD. 



3.1. Transverse Velocity of the Absorbing Material 

The facts that both of the transient features appeared suddenly, were not present in the 
immediately preceding spectra, and then persisted for the rest of the observation, imply that 
the absorbing material must have covered the entire disk of the WD within a time interval 
At of less than the exposure time of 127 s. 

With a WD radius, R^d, of 7440 ±490 km (O'Brien et al. 2001), this requires that the 
absorbing material must have traversed our line of sight to the WD at a projected speed, Vproj, 
of at least 2i?wD/At ~ 120 km s^^. Given the suddenness with which the feature appeared, 
the true velocity of the absorber, Vabs, almost certainly exceeded 120kms~^, perhaps by a 
considerable factor. 



3.2. Lower Limit on the Linear Extent of the Absorbing Structure 

The transient absorption on 1994 Oct 20 persisted until the end of the 36-min observing 
sequence (see Figs. 1 and 2). Thus the duration of the absorption event, Atabs, was at least 
1600 s. Given that the absorbing feature moved across the line of sight to the WD with a 
speed of V^bs; we conclude that the extent of the absorber in projection, AL, was at least 
Vahs Atabs, i-e., AL > 1.9 x 10^ km. Thus, the absorbing feature had a transverse dimension, 
-^CME, that exceeded ~O.28i?0. 

The dK2 star in V471 Tau has a radius that is slightly less than that of the Sun: 
i?(jK — 0.96 i?0 (O'Brien et al. 2001). Thus, the above lower limit on Lcme in terms of the 
dK star radius is of order 0.29 i?dK- 

3.3. Probability of Detecting an Absorber in Transit 

Only those CME's that emerge from the K dwarf in V471 Tau and pass in front of the 
WD will create absorption features in its spectrum. 

Moreover, the two transient features that we did detect during the 6.8 hours of observing 
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were seen to move nearly transversely across the WD line of sight in the orbital reference 
frame. 

The question we need to address is this: under what launch conditions will a CME from 
the K dwarf satisfy the above conditions? 

To answer this question, we must consider the orbital mechanics of an object that moves 
in the gravitational field of a rotating binary with mass parameter 11 — 1712/ {1111 + 1712), where 
conventionally mi is taken to refer to the more massive component. Thus for the case of 
V471 Tau we take mi to refer to the mass of the dK star and 1712 to refer to the mass 
of the WD. O'Brien et al. (2001) have used the same GHRS observations described here 
to determine the spectroscopic orbit of the WD component, and have derived masses of 
mi = 0.93 Mq and ma = 0.84 Mq. We thus have /i = 0.475. 

Using the notation of Gould (1957) for the restricted three-body problem, we adopt a 
three-dimensional {x, y, z) frame of reference rotating with the orbit, with the x-axis chosen 
to lie along the line of centers and directed from the dK primary to the WD secondary (see 
Fig. 3). The sense of rotation is denoted by the arrows close to the Q symbol in Fig. 3. With 
the origin at the center of gravity, the primary and secondary stars are fixed at Xi = —fi and 
X2 = {1 — 11), respectively. For simplicity, we restrict our attention to motions in the orbital 
{x, y) plane. 

The equations that we integrate are 

x^2y-\-x - {I- ij){x - xi)/rl - ix{x - X2)/rl , 

y = -2x + y-{l- tj)y/rl - /ly/rl 

(Gould 1957), where — Z/^ + ~ ^1,2)^- In these equations, the unit of distance is the 

separation, a, between the stars (a ~ 3.30 -R0 ~ 3.43 -Rdx), and the unit of time is such 
that Porb = Stt. The unit of velocity thus becomes 27ra/Porbj or about 320kms~^. Fol- 
lowing Gould, wc checked the numerical accuracy of each integration by evaluating Jacobi's 
constant, and found it to be invariant to within one part in 10^ after 100 time units. 

Each integration starts by launching a "test particle" radially outward from the surface 
of the K2 dwarf (with radius 0.29 in the normalized units used here). The launch point is 
chosen to lie on the equator, at a longitude 9 relative to the line of centers. Because of the 
complicated gravitational potential, the trajectory is very sensitive to both 9 and Vq, the 
initial launch velocity. For examples of the great variety of trajectories that can occur, see 
Gould (1959), especially her Fig. 8 in which fj, has the value 0.5, close to that of V471 Tau. 
If we knew how to make allowance for magnetic forces, we would probably find trajectories 
of even greater complexity. However, in the present study, we restrict attention to the 
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mechanics that are controlled by gravitation and rotation alone, which we can characterize 
with high precision in V471 Tau where the stellar and orbital parameters are well known. ^ 

We plot in Fig. 3 a series of trajectories for ejections with a range of Vq values and 
three different launch points on the dK surface. Shown as solid lines are launches from the 
sub- WD point (longitude 9 — 0), with Vq ranging from 1.2 to 2.0 (i.e., ~380-640 km s~^ for 
the parameters of V471 Tau). As noted, the plot is in the rotating frame of reference. The 
line of sight (LOS) to the Earth at the quadrature where the WD is receding from us is 
shown. The trajectories launched from the sub- WD point all pass close to the WD, and are 
favorable for creating absorption in the WD spectrum at this quadrature. Moreover, in all of 
these cases the velocity of the particle is almost entirely in the x-direction, i.e., transverse to 
our LOS. (Specifically, all of the solid hues plotted in Fig. 3 have Vy of less than ±0.07 units 
of velocity during the WD occultation. For the parameters of V471 Tau, this corresponds 
to less than ±22 km s~^.) This family of trajectories thus agrees quite well with the velocity 
properties actually observed for our two transient events. Moreover, such trajectories would 
not be able to produce absorption events at orbital phase 0.25, again in agreement with the 
fact that we did not observe any transients at that phase. 

It should also be noted that these trajectories all leave the system after passing in front 
of the WD, i.e., they represent mass that is lost from the K star and from the binary. 

Several optical observers, including Young, Skumanich, & Paylor (1988) and Bois, Lan- 
ning, & Mochnacki (1991), have reported that Ha emission tends to be strong when V471 Tau 
is near orbital phase 0.5 (i.e., when the sub- WD point faces the Earth). This might suggest 
that there is enhanced activity near the sub- WD point, which would favor launches from 
near this region on the star. On the other hand, the strong emission is probably related as 
well to the UV irradiation of the facing hemisphere by the hot WD. In any case, as noted 
by Young, Rottler, & Skumanich (1991), the Ha profiles are highly complex and variable, 
even at the same orbital phase, so that variable surface activity is clearly present over and 
above the constant effect of irradiation. 

It is possible to find other launch points and velocities that also lead to trajectories 
passing in front of the WD at phase 0.75. Fig. 3, for example, shows (as dotted lines) a 



^With increasing distance from the K star, magnetic forces on the CME wiU weaken rapidly (since 
B'^ oc r~^), whereas centrifugal forces increase in strength (oc r). As a result, the CME in effect becomes 
magnetically decoupled from the K star at a relatively short distance from the star. Beyond that point, the 
trajectory is expected to be dominated by gravitation and rotation. For our present purposes, we assume 
that the decoupling of the CME occurs so close to the K star that, apart from controlling the launch velocity, 
magnetic fields have no significant effects on the subsequent motion. 
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family of launches from 9 — +45°, with Vq — 1.4, 1.8, and 2.0, several of which do pass in 
front of the WD as desired. However, launches from 9 = —45°, shown as dot-dash lines in 
Fig. 3, do not pass in front of the WD as seen from Earth at phase 0.75. 

In order to evaluate the probabihty of CME detections such as those we observed, we 
searched in 9 and Vq to determine those launch conditions that give rise to occultations in 
which Vy is no larger than ±0.07. For purposes of occultation, wc allow the ejecta to have 
a finite extent corresponding to the lower limit mentioned above (~0.3i?dK)- as a result, 
occultation extends from slightly before to slightly after the instant at which the center of 
the cjccta transits the WD. Wc demand that Vy have values in the range ±0.07 at some 
point during this extended occultation. 

Results of this search are shown in Fig. 4, in which combinations of 9 and Vq that satisfy 
these criteria are plotted as F's. For sufficiently low Vq, the ejecta do not escape from the 
K dwarf at all; they simply fall back onto its surface. In the V471 Tau system, this happens 
at Vo,min = 1-06, shown by a vertical fine in Fig. 4. For Vq in excess of Vo,min, we find that 
launches satisfy our criteria over a certain range of 9 values. The limits of this range are 
outlined roughly by dashed lines: outside the dashed lines, the launches are not relevant 
in the present context. The allowed launch range is broadest (some 50-70° wide) for 
Vq ~ 1.4-1.6. For very fast Vq, the ejecta move past the WD at an increasingly large angle, 
such that it becomes difficult to satisfy the Vy criterion. As a result, only a narrow launch 
window (some 8-15° wide) is possible for Vq > 2. 

It is true that, in the rotating frame, the paths of most CME's eventually wrap around 
the binary: therefore many trajectories may in principle eventually intercept the LOS to 
the WD. However, this generally happens at such a great distance from the WD that the 
chances of occultation (as seen by us) are much smaller than in the closest approaches shown 
in Figs. 3-4, and also the gas density may well have fallen significantly at that time. 

If CME's with the optimal velocities Vq = 1.4-1.6 (i.e., ~450-500 km s~^) were ejected 
uniformly at all longitudes around the equator of the K2 dwarf, the probability that wc would 
detect absorption in the WD near phase 0.75 would be very roughly A^^/360° ~ 0.15-0.2. 
Thus, for every CME that is in principle detectable, we will miss A^iong ~ 5-7 others because 
they are launched from the "wrong" longitude. Mong will be even larger than this if Vq is 
outside the optimal range. 

We will also miss CME's that arc ejected out of the orbital plane. In our computer 
modelling, we considered only ejections from zero latitude. Because of the finite extent of 
Lome, however, an ejection at non-zero latitude A could still occult the WD if Lcme/o < 
sin A. With Lome ^ 0.3i?dK, and a ~ SAR^k, A exceeds ~5°. 
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If Lome were ~-RdK, A would be of order 17°. Thus, if the CME's are also launched 
from random latitudes, for every CME that is detectable, we may miss A^iat ~ 3.2-11 because 
they are launched from the "wrong" latitude. These are upper limits on iViat, since Coriolis 
forces eventually force even high-latitude ejections to approach the equatorial plane, so that 
occultation may eventually be possible. As a result, a range of crudely 3-10 is not implausible 
for iViat in the V471 Tau system. 

Thus, assuming randomly distributed launch points, for every ejection that we "catch" 
in absorption in V471 Tau at phase 0.75, we "miss" at least A^miss = Mong-Aiat ~ 15-70 
CME's. 

Of course, it is difficult to guarantee that CME's emerge from the K dwarf in a spherically 
symmetric manner. If the emergence of CME's is non-spherically symmetric, the above 
estimates of A/'miss will be modified. In the case of the Sun, where activity is confined to low 
latitudes, most CME's are observed to be launched within 40° of the equator (Hildner 1977). 
If this tendency were to occur in V471 Tau, the above estimates of A^miss should be reduced 
by factors of about 1.5. However, it is not obvious that the latitudinal distribution of activity 
on a rapidly rotating K star is the same as on the Sun: for example, such a star may differ 
from the Sun in having polar spots and in having differential rotation of opposite sign from 
that in the Sun (e.g., Vogt & Hatzcs 1991). On a star with polar spots, the CME's might 
even be ejected preferentially in the polar direction: in such a case, A^miss would be even 
larger than we estimated above. Indeed, Dopplcr mapping of the dK star in V471 Tau by 
Ramscycr, Hatzes, & Jablonski (1995) does show spots on the surface of the star at high 
latitudes. 



3.4. Rate of CME ejection 

Two transients were detected in V471 Tau in the course of 6.8 hours of our GHRS 
observations. This corresponds to a mean observable transient rate of ~7 per day. Correcting 
for the geometric factors (A^miss) in the previous subsection, and assuming that the transient 
rate that we observed is typical of the time-averaged rate, the total rate of CME emergence 
from the K star is, very crudely, of order 100-500 per day. 

3.5. Compcirison with Solar CME Properties 

We note that when CME's emerge from the Sun, the outward velocities have lower limits 
of about lOOkms"^ (Hildner 1977). This is close to the lower limit we have estimated on 
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Kbs in V471 Tau. 

However, the mean velocities of solar CME's at various epochs are considerably higher, 
350-470 km s~^ (Hundhausen 1997). In this regard, we note that the optimal launch veloci- 
ties of CME's in V471 Tau (as far as being able to create with high probabihty absorptions 
in the WD spectrum similar to those we observed) was found to be ~1.4-1.6 velocity units, 
or about 450-500 km s~^ (cf. Fig. 4). This overlaps with the high end of the solar CME 
range of velocities. 

As regards the sizes of CME's in the Sun, we note that the mean latitudinal extents of 
solar CME's at various epochs are found to be 42-47° (Hundhausen 1997). Very few events 
are larger in angular extent than 65° (Hildner 1977). Thus solar CME's near their source 
in most cases have linear dimensions that are comparable to 1 Rq, with very few having 
larger linear dimensions. If a scaling by stellar radius applies to CME's from the K star in 
V471 Tau, then Lome should not exceed about R^k in linear extent. This limit on CME 
dimensions does not disagree with the observational constraints discussed above. 

We note that inhomogeneity in solar CME's is relevant to the V471 Tau case. The 
internal structure of CME's is far from uniform. For example, two CME's reported by Wood 
et al. (1999) exhibited highly striated interiors, with many different strands of material 
apparently intersecting one another in complicated patterns. Each strand had a width of 
at least 10^ km, and was presumably associated with a different flux rope in the original 
magnetic configuration. If strands of similar widths are present in CME's from the K2 star in 
V471 Tau, then each individual strand could cover the WD (with diameter less than 10^ km) 
as a separate event. The presence of multiple strands could explain why the Si III absorption 
line fluctuates significantly with time once the event has begun, as shown in Fig. 2. 

In summary, it seems that in terms of speed, dimension, and internal complexity, the 
transient absorbing features in V471 Tau are not inconsistent with the properties of solar 
CME's. Below, we shall see that in terms of overall mass, the feature also resembles a solar 
CME. However, for the Sun, even near solar maximum, the rate of CME's does not exceed 
about 3 per day: near solar minimum, the number is closer to 1 per day (Hildner 1977). The 
rate of CME ejections from the active V471 Tau K dwarf evidently exceeds the solar value 
by a considerable factor, roughly 2 orders of magnitude. 

3.6. Another Measure of the Activity Level on V471 Tau 

Another measure of the level of magnetic activity on a star involves the size and/or 
number of cool spots on the surface. In the Sun, as the number of spots waxes and wanes 
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in the course of the 11-year cycle, the largest spots and faculae have areas of order lO"^ 
times the visible disk area. As a result of non-exact cancellation between spot darkenings 
and facular brightenings, the solar "constant" varies during a solar cycle with an amplitude 
of about 0.05%, or 0.0005 mag (Foukal 1990). 

In contrast to these small variations, the brightness of the K dwarf in V471 Tau 
varies by as much as 0.2 mag over timescales of tens of years (Ibanoglu 1978; Tunca et 
al. 1979; Ibanoglu et al. 1994). If these variations can be ascribed to waxing and waning of 
spots/faculae analogous to the solar activity cycle, then the areal dimensions of spots/faculae 
on the K dwarf in V471 Tau must exceed those in the Sun by factors of 200-400. In fact, ra- 
diometric calculations presented by O'Brien et al. (2001), along with direct Doppler imaging 
by Ramseyer et al. (1995), confirm that about 25% of the surface area of the star is covered 
by starspots, some 250 times that of the Sun. 



3.7. Dynamo Activity in V471 Tau 

We note that our estimates of the excess of the CME rates above solar values (100-500) 
and our estimates of the excess of spottcdncss above solar values (~250) are reasonably 
similar. 

Is such a high level of activity possible in the context of dynamo activity? To answer 
this, we recall that the angular velocity fl^K of the K2 dwarf in V471 Tau is some 50 times 
larger than the angular velocity of the Sun. Presumably the level of magnetic activity is 
determined by the maximum magnetic field strength on the surface of the star. Now, in a 
turbulent dynamo, the amplitude of the maximum field depends on the mechanism that gives 
rise to non-linear limiting. With one particular choice of limiting mechanism, the maximum 
field strength is predicted to scale as 17^/^ (Kippenhahn 1973). This prediction suggests that 
magnetic fields in the K2 dwarf in V471 Tau may exceed the solar values by a factor of ~350. 

We suggest that the excesses in CME rate and in spotted activity in V471 Tau are 
related to this dynamo-related excess in maximum field strength. 

What effects might arise from these strong fields? Presumably the coronal heating 
process will be highly efficient. In the Sun, enhanced magnetic field strengths, B, are cor- 
related with enhanced coronal densities, n. Because of this correlation, the Alfven speed 

Va oc Bjn^l'^ docs not necessarily vary greatly from one region of the corona to another. 
Now, when a CME is created through loss of magnetic equilibrium, the velocity of ejection 
may be related to the local Alfven speed. In such a scenario, the speeds of CME's may ex- 
hibit a well-defined mean in spite of large variations of B value from one part of the corona 
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to another. Analogously, the presence of strong fields in V471 Tau does not exclude the 
possibility that Va_ is comparable to solar Va- This may explain why the optimal ejection 
speed in V471 Tau overlaps with the mean range of solar CME speeds. 

4. Estimating the Mass-Loss Rate 

In this section, we use the properties of the transient absorption features to set limits 
on the CME mass-loss rate from the K2 dwarf. 



4.1. Column Depth of the Absorbing Feature 

We wish to estimate A^si, the number of silicon nuclei per cm^ along the line of sight 
during a transient absorption. To be specific, we consider the equivalent width of the tran- 
sient absorption in the Si III 1206 A line at the time of maximum absorption: Wa ~ 0.3 A 
(see Fig. 2). To convert this to a column density, we assume that conditions are optically 
thin; this seems to be a reasonable approximation in view of the relative weakness of the 
transient absorption. On the linear portion of the curve of growth, Wa is related to the 
column density of ground-state Si III ions, denoted A'siin, as follows (Spitzer 1968): 

Wa ^ 8.85 xlO-^nifNsiui. 

Here, Wa is in A, = 0.1206 is the wavelength in microns, and / is the oscillator strength 
of the 1206 A fine. According to Allen (1963), / = 1.9. Inserting Wa = 0.3 A, we find 
-^Siiii — 1-2 X lO^^cm"^. The total column density of sihcon nuclei, A^si, exceeds A^siiii by a 
factor l/'0(SiIII), where ■0(SiIII) is the fractional abundance of silicon in the Si III ground 
state. The numerical value of '0(SiIII) is related to the temperature of the source region 
from which the CME material emerged. 

Assuming a cosmic silicon abundance of H/Si ^ 3 x 10^ by number (Allen 1963), 
we find that the column density of hydrogen in the transient absorption is Ah ~ 4 x 
lOiV'0(Siin)cm-2. 

4.2. Number Density in the Absorbing Material 

Denoting the linear extent of the absorbing material Lqme, we recall from §3.2 that its 
value is at least ~O.3i?0 ~ 2 x 10^° cm. The volumetric number density of hydrogens, nn, 
in the absorbing material is thus given by nn ~ 2 x 10^ (O.3it!0/LcME)/'0(SiIII) cm~^. 



-14- 



We shall find below (see §4.5) that it is possible to set a lower hmit of ~10 ^ on the 
value of V'( Si III). If LcME ~ l-R©, we can therefore obtain an approximate upper limit on 
the number density of the absorbing material of nn ~ 6 x 10^^ cm~^. 

Is this limit reasonable for a detached binary system? To answer that, we note that in 
semi-detached binaries (where, due to Roche- lobe ovcrfiow, denser gas streams are expected), 
gas densities are found to be in the range 10^^-10^^ cm~^ (Batten 1970). These are, as 
expected, larger than the upper limit we obtain here for the detached V471 Tau system. 

4.3. Mass of the Transient Absorber 

With riH protons per cm^, and hnear extent of -Lome, the total mass of the 1994 Oct 20 
transient absorber, Mcme, was of order nn -^cme' where mn is the proton mass. Using the 
above estimates, we find Mcme — 3 x lO"*^^ (Lcme/0.3 i?0)^/-?/'(Si III) g. If we again assume 
LcME ~ 1 -R© (analogous to the largest solar CME's), and consider that ■(/'(Si III) < 1, then 
the CME mass for the event seen in V471 Tau must have exceeded 3 x 10^^ g. 

Is this mass reasonable for a CME? To answer this, we note that in the case of the solar 
corona, the mass distribution of CME's has an exponential distribution with typical scale 
10^^ g (Jackson & Howard 1993). Thus, once again, it is noteworthy that a CME parameter 
in the Sun has a comparable value in V471 Tau. The limiting mass of a CME from a stellar 
corona is presumably determined by the ability of the plasma and magnetic field of the 
corona to store energy only up to some hmiting value (MuUan & Mathioudakis 2000). In 
V471 Tau, the levels of magnetic activity are larger than solar, but it is less clear what the 
plasma properties might be. Thus, it is not clear whether or not the atmosphere in V471 Tau 
is in principle capable of storing larger energies than in the solar case. 

4.4. CME Mass-Loss Rate in V471 Tau 

Combining our estimates of the numbers of ejections per day (~100-500) from the K 
dwarf with Mcme from the preceding subsection, we find the CME mass-loss rate in V471 Tau 
to be of order Mcme ^ (3-15) x 10^^ (LcME/O.3i?0)V'0(SiIII) gs-^ 

How does this compare to the rate at which the Sun loses mass in the form of CME's? 
According to Jackson & Howard (1993), the maximally active Sun ejects CME mass at a rate 
of 2.7 X 10^^ gs~^. Thus the CME mass- loss rate in V471 Tau is comparable to or somewhat 
larger than that in the Sun if most of the Si is in the form of Si III, or could greatly exceed 
the solar rate if ■0(SiIII) is low. 
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We thus see that, although all of the input numbers are uncertain, the principal unknown 
factor in determining the actual CME mass-loss rate in V471 Tau is the numerical value of 
-0(81111). We now turn to a discussion of this value. 

4.5. Fractional Ionization of Silicon 

The value of ■0(SiIII) depends on the electron temperature, Te, at the place of CME 
origin. In an expanding wind, ionization states are "frozen in" within a short radial distance 

of the stellar surface. "Frozen-in" conditions set in at the radial distance where the expansion 
time of the wind becomes shorter than the ionization time-scale (Owocki 1983). Beyond that 
radial distance, even if the gas cools as the CME expands during its outward motion, the 
ionization states will not come into local equilibrium. 

We refer to ionization equilibria calculations under chromosphcric/coronal conditions 
(e.g., Arnaud & Rothenfiug 1985, their Table IV) to determine the fractional abundance of 
Si 111 at various temperatures. (At our level of approximation, we ignore the small fraction, 
generally <25%, of Si III ions that arc in excited states.) If the material in the transient 
absorber started in the corona of the K star, where electron temperatures are at least 10^ K, 
then 0(81 III) would be miniscule: Arnaud & Rothenflug's (1985) table does not extend 
beyond logTg = 5.6 for Si 111, but if we perform an (admittedly gross) extrapolation from 
the last three entries in their table, we find that at — 10^ K, ^(Si III) may be of order lO"^*^. 
If the material started out in the low transition region between corona and chromosphere, at 
(say) 140,000 K, then ^/'(Silll) ^ 10"^ The highest abundance of Si 111 occurs at electron 
temperatures of ~30,000 K, where z/'(SiIII) ^ 1. Finally, if the material started off in the 
chromosphere (at ~10,000 K) then V'(Silll) « 10"^. 

For ■0(SiIII) ~ 1, and Lcme — l-R©, the CME mass-loss rate from the results of 
the previous subsection would be (3-15) x 10^^ gs^-*^ or (5-25) x 10~^^MQyr~^. To get 
fairly reasonable agreement with the estimated lower limit of >10~^^MQyr~^ from lUE 
observations of Mg II absorption (Mullan et al. 1989), we would need -0(81111) ~ 0.005- 
0.025. These limits on t/'(SiIII) correspond to an origin either in the transition region (at 
temperatures above ~100,000 K), or conceivably in the chromosphere (at temperatures below 
~14,000 K). 

On the other hand, if the CME material in V471 Tau emerged from a truly coronal 
structure with a very low -0(81 III), then the mean mass-loss rate in CME's could be several 
orders of magnitude larger, possibly even approaching the early dynamical estimates of mass- 
loss rates as high as lO^^Mgyr"^ based on changes in the orbital period (see §1 above). 
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However, as mentioned above, those dynamical estimates failed to take into account the 
possibility of a third body in the system, or of modulations of the orbital period by the 
magnetic activity itself. 

Here we would like to point out another argument, based on accretion, which suggests 
that a mass loss rate of order 10~^ M© yr~^ is probably too large for V471 Tau. The argument 
proceeds as follows. The rate at which material is being accreted onto the surface of the WD 
has been found to be less than 10~^^ Mq yr~^ (Sion et al. 1998). As argued by Mullan et al. 
(1991) and by Sion et al. (1998), this accretion rate is several orders of magnitude smaller 
than what should be occurring if purely hydrodynamical processes were in control of the 
flow of the K2 stellar wind around the WD. This suggests that a mechanism of some kind 
(most probably a "magnetic propeller") is providing a highly efficient shield to protect the 
surface of the WD from the passing wind. The properties of such a propeller allow us to set 
a limit on the mass-loss rate. 

To determine this limit, we refer to the discussion of Mullan et al. (1991). If the mass- 
loss rate from the K2 dwarf were to be as large as 10~^ Mq yr~^, the ram pressure of this wind 
would be so large that the "magnetospheric radius," Rm, would be pushed in almost to the 
surface of the WD. Specifically, inserting a mass-loss rate of 10"^ Mq yr~^ into cqs. (13)-(18) 

1 /3 

of Mullan et al., we find that Rm would have a numerical value of only 5 x 10^ Bq cm, where 
Bq is the polar magnetic field strength in MG. Since the field strength is now known to be 
of order 0.35 MG from Zeeman splitting of the photospheric Si III lines (Sion ct al. 1998), 
Rm does not exceed 4 x 10^ cm. (The numerical value is quite insensitive to Bq because of 
the small exponent: even if Bq were to uncertain by a factor of 10, Rm would be uncertain 
only by a factor of ~2.) The above estimate of Rm is to be compared with the co-rotation 
radius, Rc — 9x 10^ cm (see eq. (20a) in Mullan et al. 1991). (The value of Rc depends only 
on the WD mass and its rotation period, both of which are rehably known.) We see that, 
with the above choice of mass- loss rate, Rm is smaller than R^. The key point can now be 
appreciated: in a situation where Rm < Rc, the propeller mechanism cannot operate. 

Since, in fact, the propeller appears to be operating with high efficiency (Sion et al. 
1998), we conclude that a mass- loss rate as large as lO^^Mgyr"^ from the K dwarf can be 
ruled out. In order to preserve propeller operation in the presence of a field of 0.35 MG, 
we find that the mass- loss rate from the K dwarf must not exceed about lO~^M0yr~^. 
This provides us with a lower limit on ■?/'(SiIII) of ~(0.5-2.5) x 10~^. Using Arnaud & 
Rothenfiug's (1985) tables, we find that the above lower limits on ■?/'(SiIII) constrain the 
electron temperature in the place of GME origin to be no more than ~230, 000-280, 000 K. 

Finally, we note that the estimates of mass-loss rate which we present here refer only 
to the component of mass loss associated with the transient events. It is possible that the 
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K2 dwarf also loses mass from its corona in the form of a steady thermally driven wind. 
But the present data have insufficient sensitivity to allow us to study the properties of that 
component directly. The upper limit on M that we derived above from propeller arguments 
refers to the total mass loss rate, i.e., a combination of steady wind plus CME's. 

5. Conclusions 

The V471 Tau system (dK2+DAl.5) provides an extremely favorable opportunity to 

detect discrete episodes of mass loss from a magnetically active star, because of the presence 
of a strongly UV-emitting, nearly point-like close companion that shines through the stellar 
wind. Magnetic activity in the K2 dwarf is kept at an unusually high level because the 
angular rotation is tidally locked to the orbital value, some 50 times solar. By extension 
from the case of the solar corona, we expect that the morphology of the corona in the K2 
dwarf in V471 Tau is dominated by magnetic structures. 

In this paper, wc report the serendipitous discovery of two transient absorption features 
of a Si III resonance line, seen in //ST/GHRS spectra. Wc interpret the transient absorption 
features in terms of material ejected from the chromosphere or corona of the K2 star and 
traversing our line of sight to the white dwarf. Our analysis indicates that the absorption 
features have properties that are comparable to those of CME's in the Sun, including size, 
velocity, complex structure, and (lower limits on) mass. The one property that shows striking 
differences from solar CME's is the rate of occurrence: we argue that the K2 dwarf in V471 
Tau ejects CME's some 100-500 times more frequently than does the Sun. 

We show that the K dwarf's mass-loss rate associated with CME's is at least (5-25) x 
lO~^^M0yr~^, but it may well be considerably higher if most of the silicon in the absorbing 
material is in ionization states other than Si III. 

In fact, the principal uncertainty in our analysis is the lack of reliable knowledge of the 
fraction of silicon, ^0 (Si III), that is in the form of Si III in the CME material; this quantity 
depends on the unknown electron temperature at the point of origin of the material. However, 
since an efficient magnetic propeller appears to be at work in V471 Tau, we can set an upper 
limit on the mass loss rate from the K2 dwarf of ~10~* Mq yr~^. 

We have recently obtained additional UV spectroscopy of V471 Tau with the Space 
Telescope Imaging Spectrograph, covering a much wider wavelength range than the GHRS 
observations reported here, and these data should provide much additional information on 
the excitation conditions in the wind and its mass-loss rate. 
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Fig. 1. — Pictorial representation of the time- resolved GHRS spectra of V471 Tau obtained 
on 1994 Oct 20, showing evidence for a coronal mass ejection from the cool component of the 
binary projected in front of the hot white dwarf companion. A portion of the UV spectrum 
centered on Lya was covered at a scale of 0.07Adiode~^. Sixteen spectra, each one having 
an exposure time of 127 s, were taken in sequence, covering a 36-min interval. The Si 111 
A1206 line suddenly appeared in absorption in the fifth subexposure and remained present 
for the remainder of the observation. Note also the broad Lya absorption wings from the 
photosphere of the white dwarf, a blue-shifted Lya emission component from the K dwarf, 
and a strong interstellar Lya absorption core. 




Fig. 2. — Equivalent width vs. time for the transient Si III 1206 A absorption feature on 
1994 Oct 20. 
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Fig. 3. — Trajectories of point masses ejected from the surface of the dK2 star in V471 Tau 
(fractional mass ratio fi = 0.4757). Only motions in the orbital {x,y) plane were considered. 
The unit of distance is the orbital separation, a ~ 3.432 RdK, between dK2 star and the 
white dwarf (WD). Solid lines: trajectories of masses launched from the sub- WD point 
(longitude 9 = 0), with launch velocities of Vq = 1.3, 1.4, 1.5, 1.6, 1.8, and 2.0 (in units of 
320 km s~^). Dotted lines: trajectories of masses launched from 9 — +45° with Vq = 1.4, 1.8, 
and 2.0. Dot-dash lines: same for 9 — —45° with Vq — 1.4, 1.6, and 2.0. The figure shows 
that launches from near the sub- WD point pass in front of the WD, as seen from Earth at 
orbital phase 0.75, with a small y-velocity, as seen in the two observed transient events. 
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Fig. 4. — Range of launch angles as a function of initial velocity, such that when the point 
mass passes in front of the white dwarf (as seen by us at phase 0.75), the radial component 
of the point mass velocity relative to the white dwarf is within the hmits ±25 km s~^ (cor- 
responding to the properties of the two observed absorption transients). Y symbols mark 
trajectories that satisfy this criterion. Launches at Vq < 1.06 (or ~340kms~^) fall back 
onto the dK star, and are marked with a vertical line labelled "no escape." 



